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ABSTRACT 

We study the electric dipole moments(EDM) of the neutron and the electron in 
the two-Higgs-doublet model, in the case that CP symmetry is violated maximally in 
the neutral Higgs sector. We take account of the Weinberg's operator 0% g = GGG as 
well as the operator O qg = qaGq for the neutron, and the Barr-Zee diagrams for the 
electron. It is found that the predicted neutron EDM could be considerably reduced 
by the destructive contribution of the two Higgs scalars to get the lower value than the 
experimental bound. As to the electron EDM, the predicted value is smaller in one 
order than the experimental one. 
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The physics of CP violation has attracted much recent attention in the light that 
the 5-factory will go on line in the near future. The central subject of the i?-factory 
is the test of the standard model(SM), in which the origin of CP violation is reduced 
to the phase in the Kobayashi-Maskawa matrix[l]. However, there has been a general 
interest in considering other approaches to CP violation since many alternate sources 
exist. The electric dipole moment (EDM) of the neutron is of central importance to 
probe these new sources of CP violation, because it is very small in SM while it can be 
larger in the other models. By begining with the papers of Weinberg[2], there has been 
considerable renewed interest in the EDM induced by CP violating neutral Higgs sector. 
Some studies[3,4,5] revealed numerically the importance of the " chromo-electric" dipole 
moment, which arises from the three-gluon operator GGG found by Weinberg [2] and 
the light quark operator qaGq introduced by Gunion and Wyler[3], in the neutral 
Higgs sector. We know the simplest extension of SM Higgs sector, namely the type II 
two-Higgs-doublet model(THDM)[6], which demonstrates explicit or spontaneous CP 
violation [7] in its neutral sector if the soft breaking term of the discrete symmetry is 
included. Many authors have already studied this model and proposed to search for 
CP violating observables directly or indirectly in the Higgs sectors[8,9,10]. 

The estimates of the maximal values of those observables are important to search 
for CP violating effect in the Higgs sector. We have already studied the Higgs potential 
to give maximal CP violation and its effect on the p parameter[ll,12]. We have found 
that the maximal CP violation is realized under the fixed values of tan (3 with two 
constraints of parameters in the Higgs potential[12]. The purpose of this paper is to 
calculate the values of EDM of the neutron and the electron in THDM with maximal CP 
violation, taking account of the contribution of the "chromo-electric" dipole moment. 
This result will give us constraints of CP violating Higgs sector by comparing with the 
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experimental bounds. 

First, we discuss the maximal CP violation in the THDM, where the Higgs potential 
with CP violating terms is written as [13]: 

v^gs = ^7i($ f i$i -hl 2 ) 2 + ^2(^2 -H 2 ) 2 

+ - |ui| 2 )($£$ 2 - m 2 ) 

+ (?'|$I$ 2 - t;> 2 | 2 + Re[h{<5>\<5> 2 - t;> 2 ) 2 ] 

+ ^ [£i _ 5al f [5i _ £al (1) 

L W 2 J [ Vi v 2 J 

where $1 and $ 2 couple with the down-quark and the up-quark sectors respectively 
and the vacuum expectation values are defined as v\ =< $5 > vac and v 2 =< $ 2 >vac- 
We do not concern ourselves here with a specific model of CP violation, but instead 
consider a general parametrization using the notation developed by Weinberg[13]. We 
take the coupling constant h in eq.(l) to be real and set 

v{v 2 = \v!V 2 \ exp(i0) , (2) 

as a phase convension. We define the neutral components of the two Higgs doublets 
using three real fields 0i, 2 , 03 and the Goldstone boson x° as follows: 

$° = ^{01 + v^ui + i{co^ X ° ~ sin/503)} , 

$° = -^{0 2 + v / 2^ + i(sm/3x° + cos/30 3 )} , (3) 

where tan/5 = v 2 /v 1 . The real fields 0i and 2 are scalar particles while 3 is pseudo- 
scalar in the limit of CP conservation. CP violation occurs via the scalar-pseudoscalar 
interference terms in the neutral Higgs mass matrix. In this basis, the mass matrix 
elements are given as 

M\ x = 20i h 2 + g'\v 2 2 + 1 \ 2 l 

\ v i\ 
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Im(/u^ 2 uf) , 



where the mass matrix is the symmetric one. Maximal CP violation is defined on a 
new basis by Georgi[14], where the Goldstone boson decouples from the $2 doublet, 
since the gauge couplings of the Higgs bosons are diagonal in this basis. The neutral 
Higgs scalars H°, H l , H 2 on this new basis are obtained by the following rotation; 
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By denoting the orthogonal matrix O that relates this basis with the mass eigenstates 
Hi, H 2 and H 3 as 

fH°\ (H 1 \ 

, (6) 





maximal CP violation is defined when 



Q2 _ q2 _ q2 
w ll _ w 12 _ ^13 



1 

3 ' 



(7) 



which was presented by Mendez and Pomarol[10]. Here, the matrix O is related with 
the orthogonal matrix U, which is defined in ref.[12] to diagonalize the mass matrix of 
eq.(4), as follows: 

cos (3 sin (3 \ 
O = I -sin/3 cos/5 U . (8) 
1/ 
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In ref.[12], the constraints of the Higgs potential parameters to get maximal CP 
violation have been studied. In this paper, two solutions yielding maximal CP violation, 
which satisfy the condition of eq.(7), for tan /3 [12] have been obtained; 

Soil : tan/5 = -L(>/3-l) = 0.51 • • • , Sol.II : tan/5 = -L(V3+1) = 1.93 • • • , (9) 
v 2 v 2 

with constraints 

(?x + ^ 2 + 2 5 - 2£ = , ffl = g 2 , = ^ , (10) 

where £ = £/|t>it>2| 2 and CP violating phase takes its maximal value as is expected. 
The condition of g\ + g 2 + 2g — 2£ = gives an important constraint for the neutral 
Higgs scalars and the charged Higgs one, since £ determines the charged Higgs mass as 
follows: 

m 2 H± = lv 2 , (11) 

where v 2 = vf + v 2 . In addition to these constraints, there are the positivity conditions 
such as[15] 

9i > 0, 52 > 0, g > -y/gig2, 9 + g' - \h\ > -y/glg*, 

£ >0, g'-\h\+l>0, l-g>-^g~ 2 . (12) 

The masses of the three neutral Higgs scalars are given as 



rn 



HI 



= 20i cos 4 p + 2g 2 sin 4 + 4(f - sin 2 /3 cos 2 /3 = 2 gi , 



™H2 = + £ + ™L = </ + £-^, (13) 

in units of v 2 , where the conditions in eq.(10) are used in the second equality for m 2 H1 . 
We notice that the four Higgs masses are given by four parameters g±, g', h and £. Since 
the parameter h is predicted to be very small in analyses using the renormalization 
group equation[16], the values of m H2 and m H3 are expected to be close to each other. 



In our numerical analyses, we take mni < rnm by fixing h < as our convention. On 
the other hand, mux is not constrained. 

Now, we can estimate the CP violating parameters lmZ 1 and lmZ 2 in THDM. 
These are the imaginary parts of the scalar meson fields normalization constants, Z^ 
which are the column vectors in the neutral Higgs scalar vector space, defined in terms 
of the tree level approximation to the two-point function as follows[13]: 

i<^ = E^%"^ B) (^ = 1 ' 2 )' (14) 

v i n =l 5 _ m Hn 

where Vi = (0°) mc . The CP violation factors ImZ^ are deduced to 

( k ) tan/3 (fc) 1 (fc ) (fe ) 

ImZ} = -w 'u^ , ImZ^ = — — -u\ >u' 3 , 15 

cos p tan p sin p 

where uf^ denotes the i— th component of the A;— th normalized eigenvector of the Higgs 
mass matrix. The values of uf^ are given by solving Higgs mass matrix M 2 of eq. (4) 
in the case of the maximal CP violation. In the maximal CP violation, we get [12] 

= (cos/5, -sin/3, 0) , 
vP^ = (— ^=sin/3, — ^=cos,3, — \=) , 

m (3) = (-L sin/?, -^cos/3, -=) . (16) 

Using these eigenvectors, we can calculate CP violating parameters ImZ; . For the 
first Higgs scalar, these are zero because the third component of the eigenvector is 
zero as seen in eq.(16), i.e., there is no scalar-pseudoscalar interference term. We have 
non-vanishing values for second and third Higgs scalars (k=2,3) as follows: 

lmZ[ 2) = -lmZ[ 3) = i(V3 T l) 2 , ImZf = -ImZ^ = ~^{V3 ± l) 2 , (17) 

the upper and lower signs correspond to the Sol. I and Sol. II of tan/? in eq.(9), respec- 
tively. We notice that these values are somewhat smaller than the Weinberg's bound[13] 
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taking the same value of tan (3, 




{ 



0.89 
0.46 




0.46 
0.89 



(18) 



where (WB) denotes the Weinberg's bounds, and the upper values and lower ones 
correspond to the Sol. I and Sol. II of tan f3, respectively. Thus, the Weinberg's bound 
does not correspond to maximal CP violation, because both ImZ[ k ^ and ImZ^ cannot 
approach to those bounds at the same time. 

Let us discuss the EDM of the neutron in the case of maximal CP violation. The 
low energy CP violating interaction is described by an effective Lagrangian Lcp, which 
is generally decomposed into the local composite operators Oj of the quark and gluon 
fields, 



Some authors pointed out [2,3] that the three gluon operator with the dimension six 
and the quark-gluon operator with the dimension five dominate EDM of the neutron 
in THDM. So, we study the effect of these two operators on the neutron EDM. Var- 
ious techniques have been developed to estimate the strong-interaction hadronic ef- 
fects[17,18,19,20]. The simplest one is the NDA approach[17], but it provides at best 
the order-of-magnitude estimates. The systematic technique has been given by Chem- 
tob[18] for the case of the operator with the higher-dimension involving the gluon fields. 
We employ his technique to get the hadronic matrix elements of our operators. 

Let us define the following operators: 



where q denotes u, d or s quark. The QCD corrected coefficients are given by the 



(19) 



O qg {x) = -^qa^q , 3g (x) = -|f 




(20) 
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two- loop calculations [2, 3] as follows: 

y/2G F m u {n) 2 / g s (/i) ■ 

C d9 = _ ^^(^) g 2 M[/( , t)tan2/3ImZ2 _ g^cot^ImZi] (jM 
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C 3 , = ^ M , ()ImZ2 (^ , (21) 

where z t = {m t /m H ) 2 and we omitt the upper- indices (A;) defined in eq.(15). The 
functions f(zt), g{zt) and h(zt) are the two-loop integral functions, which are defined 
in refs. [4,5,21]. The coefficient C sg is obtained from d g by repairing rrid with m s . 

The hadronic matrix elements of the two operators are approximated by the inter- 
mediate states with the single nucleon pole and the nucleon plus one pion. Then, the 
nucleon matrix elements are defined as[18] 

{N(P)\Oi(0)\N(P)) = AiU(P)i>y 5 U(P) , 
(7V(P')|0^(P)7r(^)) = B i U(P')T a U(P) , (22) 

where U(P) is the normalized nucleon Dirac spinors with the four momuntum P. By 
using Ai and B^i = ug, dg, sg, 3g), the neutron EDM, d^, is written as 

d l = 7TT C i A i + f (9ttNN, m n , m n ) ^ c * B i > ( 23 ) 

Zm n i i 

where /i n is the neutron anomalous magnetic moment. The function F^tvv, "^n, 
was derived by calculating the pion and nucleon loop corrections using the chiral La- 
grangian for the coupled Nn^ and is given in Appendix A of ref. [18]. Here, the dimen- 



sional regularization with the standard MS scheme is used for defining the finite parts 
of the divergent integrals. The coefficients Ai and E>i were given by the use of the large 
N c current algebra and the i] meson dominance [18]. Then, we have 
A - fn R 4(m u + m d )q 1 / i 

A i - JiQrioNN , Di — — , (24J 



with ai = -(m E o - m s )/(2m s - m u - m d ) ~ -0.28 and = ^2/3F = 0.186GeV, 
where fa is defined by 

Mg)|O,(0)|0> = f\q 2 . (25) 
The values of fi were derived by using QCD sum rules as follows[18]: 

f qg = -0.346GeV 2 , f 3g = -0.842GeV 3 , (26) 

where f qg denotes the flavor singlet coupling. 

Now, we can calculate the neutron EDM. Our input parameters are [20] 

A QCD = 0.26GeV , (m u , m d , m s ) = (5.6, 9.9, 200)MeV , p = m n , 

M = m t = 150GeV , g nNN = 13.5 , g VoNN = 0.892 . (27) 

We show in fig.l our predictions of the neutron EDM versus the ratio of the masses 
of the two Higgs scalars in the case of m H2 = 200, 400, 600GeV by using Sol. I of 
the maximal CP violation. The predicted neutron EDM in Sol. II is almost same as 
the one in Sol. I. As seen in fig.l, the mass difference Amj? = muz — Tan2 should be 
small in order to get the lower values of EDM than the experimental upper bound, 
11 x 10~ 26 e-cm[22]. We get Am H ~ 18GeV, which corresponds to \h\ ~ 0.13 , in the 
case of m H2 = 200GeV, Am H ~ 52GeV(|/i| ~ 0.73) in the case of m H2 = 400GeV, and 
Am;; ~ 140GeV(|/i| ~ 3.12) in the case of m H 2 = 600GeV. In order to show the effect 
of the relevant operators, we show the contributions of O ug , Od g + O sg and 3g to the 
neutron EDM in fig.2 in the case of m# 2 = 200GeV by using Sol. I. As seen in fig. 2, 
the effects of the O ug and 3g cancel out each other, and the Od g + O sg contribution 
dominates the neutron EDM, which has already been pointed out in ref.[20]. Since the 
contributions of O ug and 3g depend only on lmZ 2 as seen in eq.(21), there are large 
difference between Sol. I and Sol. II in these contributions. However, these contributions 
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Figure 1: The predicted neutron EDM versus the ratio of A^H3 in the cases of 
m>H2 = 200GeV(solid curve), 400GeV(dashed curve) and 600GeV(dashed-dotted curve) 
by using Soli. The horizontal dashed line denotes the experimental upper bound, 11 x 
l(T 26 e • cm. 



almost cancel out each other in the predicted total neutron EDM. On the other hand, 
the contribution of Od g + O sg in Sol. I is exactly same as the one in Sol. II. Thus, both 
Sol. I and Sol. II give the almost same values for the predicted neutron EDM. 

Now let us discuss the EDM of the electron. Barr and Zee [5] presented the two- 
loop Feynman diagrams which can lead to a large EDM of the charged lepton. Those 
diagrams involve a heavy particle, say the top quark or W boson in the loop that 
couples to an external photon line as follows: 



d, 



= - ^V2G F m e {f(z t ) tan 2 f31mZ 2 + g(z t ) cot 2 pimZ^ , (28) 

t—loop 



d, 



a 



W-loop 327r3 



\f2G F m e [Sf(z w ) + hg(z w )\ (sin 2 (3 tan 2 plmZ 2 + cos 2 pimZj , 



where z t = (m t /mH) 2 and Zw = (jnw / m ii) 2 ■ The loop-functions f(z) and g{z) are the 
same ones as those in eq.(21). Chang, Keung and Yuan[8] have given the complete set of 
two- loop diagrams in the multi-Higgs-doublet model. They calculate the effective CP 
violating HZy vertex in addition to the H^yy one induced by the unphysical charged 
Higgs and the W contribution. In this paper, we use their results[8] instead of eq.(28), 
however the structures of tan/? and ImZi in eq.(28) being unchanged. 

In the case of maximal CP violation, the loop contributions of W and the unphysical 
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Figure 2: The components of the predicted neutron EDM versus the ratio of vtiu^I^hj, 
in the cases of m^ 2 = 200GeV by using Sol. I. The thick-solid curve denotes the abso- 
lute values of the total contribution. The thin-solid curve, the dashed curve and the 
dashed-dotted curve denote the contribution from O ug , Od g + O sg and 3g operators, 
respectively. 
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Figure 3: The predicted electron EDM versus the ratio of m^/miis i n the cases 
of rriH2 = 200GeV(solid curve), 400GeV(dashed curve) and 600GeV(dashed-dotted 
curve). The experimental upper bound (—0.3 ± 0.8) x 10~ 26 e • cm is outside the figure. 



charged Higgs boson vanish because sin 2 /? tan 2 /3lmZ 2 + cos 2 filmZx is zero in both Sol. I 
and Sol. II. So, the elctron EDM is given by only the top-quark loop contribution, which 
is exactly same in both Sol. I and Sol. II. We show in fig. 3 the predicted electron EDM 
versus the ratio of the two Higgs scalars in the case of m H2 = 200, 400, 600GeV. 
Thus, the magnitude of the electron EDM does not go over the experimental bound 
(—0.3 ± 0.8) x 10~ 26 e-cm[22] even if the mass difference of the two Higgs scalars is 
considerably large. Summary is given as follows. We have studied the EDM of the 
neutron and the electron in the two Higgs doublet model, in the case of CP symmetry 
being violated maximally in the neutral Higgs sector. The maximal CP violation is 
realized under the fixed values of tan (5 with two constraints of parameters in the Higgs 
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potential. We have taken account of the Weinberg's operator 3g = GGG and the 
operator O qg = qaGq for the neutron, and Barr-Zee diagrams for the electron. It 
is found that the predicted neutron EDM could be considerably reduced due to the 
destructive contribution of the two Higgs scalars leading to the lower value than the 
experimental upper bound, while the predicted electron EDM is smaller in one order 
than the experimental bound. Since our predicted value of the neutron EDM lies 
around the present experimental bound, we expect its experimental improvement to 
reveal the new physics beyond SM. 
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